What happens to turbulent skin-friction drag reduction at high Re? 
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We address one of the capital problems in the field of turbulent skin-friction drag reduction, i.e. the 
performance of the known techniques at high values of the Reynolds number Re. We limit ourselves 
to considering open-loop techniques based on spanwisc forcing (oscillating wall, streamwise-travclling 
waves), and explore via Direct Numerical Simulations (DNS) how quickly the drag reduction and 
net energy savings decrease when the friction Reynolds number is increased. We consider turbulent 
channel flow up to Re T = 2000. Owing to the demanding computational cost of large-scale simula- 
tions, the size of the computational domain in the homogeneous directions is reduced while keeping 
under control the statistical uncertainty related to the finite time average. This allows us to carry 
out a full parametric study (at least at Re T = 1000), hence avoiding the two assumptions that the 
low-i?e optimum set of parameters remains optimal at higher Re, and that the optimum obeys an 
inner scaling. 

Although this computational strategy definitely decreases the absolute fidelity of the DNS results, 
a careful estimation of the errors allows us to interpret the obtained results by suggesting an unex- 
pected and interesting scenario where the performance of the drag-reduction technique degrade with 
Re with a markedly different rate depending on the parameters. In particular, the known optimal 
region turns out to be such at low-i?e only, since there drag reduction degrades quite fast with Re, 
in line with available results. However, other regions arc much less sensitive to Re, or insensitive 
at all. If one considers that the energy required to create the forcing presents a slightly favorable 
trend with Re, the possibility exists of net energy saving at very high Reynolds numbers. 

This interesting scenario remains speculative in nature, owing to the spatial truncation implied 
by the limited domain size. However, a few full-scale DNS for the traveling waves at Re T = 400 have 
been carried out, and their results fully support the suggested scenario, which, though appealing, 
would force us to rethink our current understanding of how these drag reduction techniques work 
and behave at high Re. 
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I. INTRODUCTION 

Decreasing the aerodynamic drag in the turbulent regime is a technological objective that is prompting 
the study and development of several techniques aimed at skin-friction drag reduction, being the increase of 
friction drag above the laminar value one of the fundamental manifestations of turbulence in the simplest 
wall-bounded flow. The techniques under development possess a wide range of effectiveness, ranging from 
passive techniques (a classic example being ribletsi^) that yield up 8-10% drag reduction in well controlled, 
low-Reynolds laboratory conditions, to active techniques based on linear control theory 2 ^, that promise much 
better performance, especially in terms of net energy savings, but are characterized by high complexity 
and have been studied through high-fidelity Direct Numerical Simulations (DNS). Midway between these 
two extrema, open-loop active techniques compromise between energy gain and energy expenditure to yield 
sizable net energy savings with the advantage of moderate complexity, as they do not require the deployment 
of micron-sized sensors and only require relatively large-scale actuators. One recently proposed strategy, 
that looks promising in terms of net energy saving potential, is the streamwise-traveling waves concept 2 ^, 
that has been shown to provide at least 20% net energy savings and has been already given experimental 
verification^. The long-term perspectives of these drag reduction techniques are illustrated in a recent 
volume devoted to the subject 2 ^. 
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Besides the obvious need for further technological development of sensors and actuators, essentially all the 
active techniques lack a thorough evaluation in flows characterized by high values of the Reynolds number. 
Unfortunately, however, most of the anticipated applications are at high Re. In this paper we critically 
address the effect of varying the Reynolds number on wall-based open-loop control techniques by focusing 
on one particular control strategy: the streamwise-traveling waves of spanwise wall velocity 



where W(x, t) is the velocity forcing at the wall, of amplitude A. This forcing law contains the two cases of 
the spanwise-oscillating wall, for A x — ► oo, and the steady wave when u = 0. 

It is the present Authors' opinion that this is becoming the most fundamental question that should be 
answered in this research field. The reason why few data exist at high Re are quite obvious, and can be 
traced back to the enormous increase of the computational cost of DNS for increasing Re, paralleled by 
a shrinking of the physical size of the required sensors and actuators to be employed in an experiment. 
Using wall-based strategies to affect the turbulence wall cycle, that was shown to be autonomous^!, does 
not guarantee that performances scale up at high Re, and in fact the few available data, covering less than 
one decade in Re, witness a significant decrease of, for example, maximum drag reduction. It remains to 
be understood whether such a decrease is a low-i?e effect, so that performances are deemed to stabilize 
as soon as Re becomes large enough, or is expected to continue indefinitely, so that the level of zero drag 
reduction is reached at values of Re which are still smaller than those typical of applications. Figure [T] 
graphically illustrates the two scenarios and shows that both are compatible with the available low-i?e 
numerical and laboratory experiments, so that it can not be inferred which one is more likely to occur. 
The only information available is the result of the study by Iwamoto et ali^, who used both analytical 
developments and DNS experiments to show that a virtual drag reduction technique capable of removing 
near-wall turbulent fluctuations in a turbulent plane channel flow for y + < 10 still yields 35% drag reduction 
at the relatively large value of Re T = 10 5 . 

We are unaware of studies where the spanwise-traveling waves by Du & Karniadakis* have been studied 
as a function of the Reynolds number. For the present class of techniques, the available information on the 
high-i?e behavior is as follows. J.I.Choi et al^ used turbulent channel flow simulations at Re T = 100, 200 
and 400 to suggest that, at T + = 100 and A + = 10, the oscillating wall produces a drag reduction of 
41.1%, 29.9% and 22.4% at the three Reynolds numbers respectively. They fitted their results with those 
available in the literature to infer that the drag reduction is predictable by a quadratic function of a scaling 
parameter V c that is proportional to i?e~ 0,2 . The finding is coherent with the results reported above that 
show a Re~ 0A dependency. Ricco & Quadrio^S reported that the spanwise-oscillating wall at T + = 125 and 
A+ = 12 yields 32.5% drag reduction at Re T = 200 and 28.1% at Re T = 400. Compared to the results of 
J.I.Choi et al^, a slower decay against Re is obtained, of only Re~ 02 , probably owing to the differences in 
the forcing parameters and in the numerical discretization, being the resolution higher in the latter case. 
Quadrio, Ricco & Viotti 2 ^ determined the best-performing traveling wave at Re T = 200 and forcing intensity 
of A + = 12, and verified that at Re T = 400 maximum drag reduction decreases from 48% to 42%, again 
supposing a power law dependency, R would decrease as Re~ 019 . Touber & Leschziner— present DNS data 
at Re T = 200, 500 as well as LES and one DNS datapoint at Re T = 1000 to propose that drag reduction 
values obtained under similar values of a scaling parameter S decrease roughly as Re~ ' 2 , thus confirming 
the findings of J.I.Choi et al^. Lockerby et ali^ presented a linear stability study to investigate the influence 
of a spanwise generalized Stokes' layer on the growth of near- wall turbulent streaks. The streak amplification 
is shown to be strictly related to the drag reduction predicted by numerical experiments, suggesting the 
importance of basic linear mechanisms on the drag reduction. Their linearized model allows predictions 
at higher Reynolds number, and the effect of the Reynolds number is found to be much milder than that 
obtained in simulations. In a recent paper by Moarrcf & Jovanovio2£, a model-based approach that employs 
linearized Navier-Stokes equations in conjunction with turbulence modeling is developed. The i?e-effect up 
to Re T = 934 is addressed for the spanwise-oscillating wall: the optimal frequency of oscillation is found 
to increase with Re, and the drag reduction decreases as i?e~ 015 , i.e. again at a slightly slower yet still 
comparable rate than numerical experiments. If exception is made for the last two works, that span a 
parameter range but are limited by their inherent uncertainty, all the others tracked only the neighborhood 
of the forcing parameters yielding maximum drag reduction at low Re, for the oscillating wall case, on 
the assumption that these optimum forcing conditions scale in wall units. We are unaware of laboratory 
experiments that have been able to significantly vary the value of Re. Although it is true that experiments 
are typically carried out at values of Re which are somewhat higher than those typical of DNS (typically 
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FIG. 1. Variation with Re of maximum drag reduction R ma x for spanwise forcing as educed from several studies. 
Black symbols indicate numerical (DNS) studies, and open symbols are for experimental studies. We explicitly note 
that the forcing amplitude is not always comparable among different datasets. Circles: oscillating waH^^ ' 18 ' 24 ' 27 ' 29 ^ ; 

1 28 S 9 3 2^ 

triangles: streamwise-travelling waves-^; squares: spanwise-travelling waves-—; triangles: Lorentz force-^. The 
gray area marks the "unknown zone" where drag reduction could drop according to one of the two lines drawn in 
the sketch, i.e. falling rapidly to zero (blue line, scenario 1) or gradually decreasing the drop rate (red line, scenario 
2). 



Re T = 600), experimental uncertainty, and above all the issue of the spatial transient that makes a control 
surface only partially effective precludes putting together numerical and experimental data to precisely assess 
trends of DR reduction with Re based on mixed datasets. 

In this paper, we use DNS to obtain information on the performance of the spanwisc-oscillating wall 
and the streamwise-travcling waves up to Re T = 2000, using DNS. The aim is to be able to carry out a 
parametric study at higher Re, so as to avoid the double hypothesis implicitly assumed in the previous works. 
In fact it has traditionally been assumed that at low-i?e optimal parameters determined at low Re remain 
optimal at higher Re, assuming they scale in wall units. To keep the computational study manageable, 
our strategy is one of changing the size of the computational domain while carefully controlling the effect 
of this essential discretization parameter on the reliability of the obtained drag reduction information. On 
one hand, random and systematic errors may be introduced and are rigorously estimated. On the other, 
the simplified framework of a narrower domain can give useful suggestive results on the behavior of drag 
reduction and on the importance of the truncated large turbulent scales on it. 



The structure of the paper is as follows. Firstly the numerical strategy and simulation parameters are 
described in detail, focusing on the eventual convenience of using a reduced domain, on the strategy we 
adopted to estimate random errors on the mean skin-friction. The parameters needed to evaluate the 
control are described as well before estimating the influences of a smaller domain on the skin- friction of the 
uncontrolled and controlled channel-flow. In the third section the drag reduction ratio for traveling- waves up 
to Re T = 2000 is discussed and the so-called .Re-effect quantified. In the fourth section interesting insights 
on the power consumption and power budget are given, followed by concluding remarks. 
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TABLE I. Computational domain size and friction Reynolds number Re T - 



II. METHOD 

A. Code and simulation parameters 

The analysis is carried out on a newly created DNS datasct. The computer code used for the calculations, 
its parallel algorithm and the architecture of the computing system have been described elsewhere^. The 
code is a mixed-discretization parallel solver of the incompressible Navier-Stokes equations, based on Fourier 
expansions in the homogeneous directions and high-order explicit compact finite-difference schemes in the 
wall-normal direction. Most of the calculations have been run on a computing system available in dedicated 
mode at the Universita di Salerno. If exception is made for the domain size, discussed below, there is 
nothing specific about the present set of calculations. The governing equations are integrated forward in 
time, starting from an initial condition of fully developed uncontrolled channel flow, generated specifically 
for each Reynolds number. The flow rate is kept constant and the bulk velocity is Ub = 2/3U p , where U p is 
the centerline velocity of a Poiseuille flow with the same flow rate. Parametric DNS studies have been carried 
out at Re p = U p h/v = 4760 and 23500, corresponding to Re T = u r h/v = 200 and 1000 respectively, based 
on the friction velocity u T of the unmanipulated flow. A more limited dataset is produced at Re p = 73000, 
corresponding to Re T = 2000, the highest Reynolds number reached during this investigation. 

The spatial resolution in the stream- and spanwisc directions is Aa; + = 10 and Az + = 5, the wall-normal 
resolution varies from Ay^ = 1 at the wall to Ay+ a:t , = 6.9 at the centerline, with the number of grid 
points in this directions set at Re T /2. The resolution has been kept constant in wall units while varying the 
Reynolds number. In all of the traveling-wave calculations at Re T = 200 and Re T = 1000 the streamwise 
length was set to L+ = 1250 wall units, while the spanwise length was Lf = 625. The domain size has 
been increased to L+ = 2700 and Lf = 1350 at Re T = 2000 to improve the reliability of the results. We 
note that such domains are of reduced sizes compared to the computational domains that are customarily 
employed but are still rather large when compared to the Minimal Flow Uniti£, the smallest domain that 
can self-sustain turbulence. The domain used for the oscillating wall (Rcfl in tabic U ) is slighly smaller, 
being L+ = 2L+ w 1000. As figure [3] shows, the domains used for all the traveling-wave cases are only three 
times smaller in width and length, compared to a typical domain used in channel flow simulations, while 
they are still 20 times bigger than the minimal domain. The domain sizes used for most of the simulations 
at each Reynolds number are shown in bold in tabic HI together with other domain dimensions that have 
been simulated to run the size-sensitivity tests discussed below. 

These box sizes are the result of an accurate a posteriori study to find the best compromise among 
accuracy, reliability and computational efficiency. Nonetheless, the following a priori criteria identify a 
bound for the domain extent required to represent the drag reduction phenomenon reasonably well. 

First of all, the domain should be large enough to avoid any obviously nonphysical narrow-channel effect, 
such as one-sided turbulence^. This basic requirement is always amply verified in our simulations. Next, 
all the turbulent structures that affect the wall shear and are thus potentially involved in the drag reduction 
should be present, but this is a more stringent requirement as well as more difficult to be clearly stated. 
Although the periodic boundary conditions imply an in(dc)finitc spatial domain, choosing the fundamental 
wavenumbers in the horizontal directions sets an upper limit for the structures that fit into the domain. 
In the near-wall region, the size of the turbulent structures increases with the distance from the wall. By 
increasing the distance form the wall, the box width or length may become too small, so that the size of 
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FIG. 2. Size of a regular domain compared to the reduced domain used in his work and the Minimal Flow Unit 
presented by Jimenez and Moin^. 

the structures is constrained loosing its dependency on the wall distance. Based on the aspect-ratio of the 
most common near-wall turbulent structures, the turbulence will be well-resolved up to a distance yd from 
the wall if the spanwise length L+ « . The final requirement on L z is found by considering that the 

logarithmic region has to be well resolved to obtain good prediction of the wall shear—. 

The streamwise length of the computational box has been dictated by the fundamental wavelength of 
the wall forcing, or by the will to keep an x — z aspect ratio of 2:1. In any case, the minimal fimrbai for a 
good description of the logarithmic layer L+ ps h + has been respected. 

A real benefit is obtained by the use of reduced domain when the averaging time is adjusted, so that a 
slight uncertainty on the wall shear prediction is still allowed. This will keep under control the simulation 
time and hence the computational cost. After various attempts, simulations were run for a time equal to 
TsimUp/h = 1500 at Re T = 200, in viscous timescale Tj im = 12500. At Re T = 1000 simulation time has been 
kept constant in viscous unit, leading to a simulation time T S i m U p /h = 400. As a result, the uncertainty 
level on wall shear is increased, because both the actual averaging area and integration time in outer units 
arc reduced. At Re T = 2000 a time of T S i, n U p /h — 400 has been used, resulting in T+ m = 24000, in order 
to increase the reliability of the results. 

The control parameters are chosen such as to to compare present results to those presented by Quadrio 
et a.&. In particular, A + = 12 is chosen, although this is not the optimal choice to study the net power. 
For clarity sake we report in figure |3] the drag reduction ratio results by Quadrio et al2£ accompanied by 
five dashed lines marking the regions of the parameter space investigated in the present work. Two scans 
at constant wavenumber have been made, one for the oscillating wall case at k x = (line 1) and one at 
k x = 0.005, corresponding to A+ = 1250 (line 2). Here, more full channel DNS data are available and 
the maximum DR is achieved at low-i?e. Two more sweeps at constant frequency uj + = —0.2 (line 3) and 
uj + = 0.012 (line 4) have been done to study the effect of a change in wavelength. Several waves on the 
locus of maximum DR, at constant wavespeed c + = 2.5, are also considered (line 5). A wide region of the 
parameter space is then spanned. The maximum frequency and wavenumber analyzed are equal, in wall 
units, to those of Quadrio et al^. 



B. Computational cost 

The existence of a long initial transient during the onset of statistically-steady drag-reduction is a peculiar 
feature to the present problem that makes the run of parametric studies particularly expensive. Its duration 
does not depend on the domain size and for a full DNS it usually takes a significant part of the computational 
effort only to get rid of the initial transient. This is crucial when trying to achieve a full parametric study 
of an open-loop flow-control technique. 

The simulation of smaller domains clearly reduces the time needed to advance the simulation by one time 
step and thus the time spent in resolving the transient. The per-timestep computational cost due to the 
pseudo-spectral convolutions can be expressed, following Jimenez^, as Q ~ N y N x N z \og(N x N z ), where N x , 
N y and N z are the number of modes/collocation points in the three directions. If we consider a domain 
aspect ratio of L x /L z = 2/1 and a resolution Ax + = 2Az + = 10, the number of modes in the homogeneous 
direction is N x = N z = L+/Aa- + . It follows that the time required to advance a timestep is more than 16 
times shorter if a domain box that is four times shorter and narrower than a regular domain is used. 
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FIG. 3. Map of the drag reduction ratio 100R versus the wavenumber n x and frequency cu of the traveling wave at 
A + = 12 and Re T = 200, after Quadrio et The dashed lines represent the portions of the parameters space 

investigated in this work. Dark red: constant k x simulations; black: constant cj simulations; white: constant c 
simulations. 



If on one hand the wall-clock time needed to reach a statistically-steady flow is shortened, on the other 
hand with a restricted domain the averaging time needs to be larger than 10 wash-outs to still yield "good" 
statistics, i.e. statistics with the same uncertainty level (of course here we are neglecting any possible 
narrow-channel effect). The standard deviation of the mean wall-shear can be expressed as: 

a T — cr T /y/T sim (1) 

where a T is the standard deviation of the spatially-averaged wall-shear signal that is proportional to 
(A + )~ 1 ' 2 , where A — L x x L z is the area of the domain. This means that a 16 times smaller domain 
would require a 16 times longer average to keep the same uncertainty level. It is then clear that an advan- 
tage can be achieved with a smaller domain mainly if the domain size is decreased while keeping a finite 
averaging time. We allowed then for a certain value of uncertainty on the measured skin-friction that clearly 
needs to be estimated in terms of confidence interval in order to have a meaningful result. 



C. Random uncertainty evaluation 

A strategy is needed to estimate the uncertainty on the measured drag reduction ratio R that involves 
only the time-resolved data of the space-averaged wall-shear stress r. When a finite averaging-time is used 
in combination with a domain of reduced size, <r r increases while T s - lm is finite, so that the uncertainty on 
the mean due to the intrinsic randomness of the process may not be negligible. 

The method used in this work is based on the assumption that statistical data are single realizations, 
equally spaced in time, of a continuous stationary (ergodic) random processes. Strong stationarity specifies 
that all statistical moments are time-invariant. Implicit in this is also the assumption that the signals are 
free from any initial transient. Under these hypotheses, the standard deviation of the mean oy is related to 
the autocorrelation function p(At) as follows: 
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FIG. 4. Values of the friction coefficient Cj versus the size L X L Z of the computational domain expressed in wall 
units. Horizontal lines are the values of Cj as educed from the Dean correlation. 
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FIG. 5. Values of the drag reduction 100i? achieved with the oscillating wall, at Re T = 200, T + = 125 and A + = 12, 
versus the size of the computational domain expressed in wall units. The horizontal line marks the value computed 
by DNS in Quadrio & Ricco^i. 
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where p(At) = E["r(i)T(£ + At)]/al and E[-] is the expected value operator. When T s ; m 3> At max , then the 



expression reduces to a- = 2(cr^/T sim )T, where T — Jq' 1 " 1 p(At)dAt is the integral timescale of the process 
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The timescale accounts for the need to discard data points in the time-resolved history that can not be 
considered uncorrelated. The variance of the wall shear and the autocovariance involved in the definition of 
autocorrelation are themselves unknown and can not be measured directly, therefore their best estimators, 
the sample standard deviation squared and the sample covariance q(At), are used instead. 

The standard uncertainty of the drag reduction ratio sr is computed by propagating via RSS the sample 
standard deviation of the mean skin-friction of both the uncontrolled s T o and the controlled Sr case. A 
confidence interval for R can be obtained thanks to the central limit theorem, hence: 

R - 7«R < Mr < R + 7 S R 

where 7 = z a / 2 is the standardized confidence interval of a normalized Gaussian PDF and depends on the 
desired confidence level 1 — a. In this work the common confidence level of 0.95 has been chosen. 



D. Parameters for control evaluation 



The control performances are evaluated through the three dimensionless parameters introduced by Kasagi 
et alJ^, (i?, Pi n and S), that describe the total energy budget. The time- averaged pumping power per 
channel wall area is defined as: 



T ■ —T- 



T sim 

<t x > dt 

Ti 



where t x denotes the streamwise component of the wall-shear stress and < • > is the wall space-average 
operator. The reduction of the pumping power due to the control is then 

R=^ 

where P and Po are the pumping powers required to drive the controlled and reference flow respectively. 
It is important to notice that in a constant flow-rate simulation the reduction of the pumping power is 
identically equal to the reduction of the skin-friction drag. The power introduced into the system to apply 
the control is the input power per unit area Pin-, defied as the power needed to drive the wall against the 
spanwise wall-shear: 

p in = E 7^- ^ / TS "" < Wt z > dt (2) 



Pq (?si m — T{ 



Ti 



where r z is the spanwise component of the wall-shear stress. The budget between benefits and costs of the 
control can be then easily defined as: 

S = R Pin 



E. Narrow-channel effects 



The use of smaller domains brings about a potential bias in the flow statistic and, in particular, in 
the mean skin-friction coefficients. It turns out that the skin-friction is increasingly underestimated for 
decreasing domain sizes. The effect is evident in figure 21 the relative error between the computed C/ and 
that predicted by the Dean's^ correlation is shown against the domain size in outer units. A deterministic 
behavior is found and in this range of Re an empiric power-law holds between the relative error and the 
domain area. At the larger domain sizes used at Re T = 200 the trend of the relative error is shaded by the 
statistical error and by the known slight inaccuracy of the Dean's correlation at low Re. The domain size 
at Re T 1000 and 2000 has always been chosen to have a relative error of less than 10%. 

The drag reduction ratio is itself affected by a change in the domain size. Even though at Re T = 200 the 
mean skin-friction coefficient of the uncontrolled channel is almost insensitive to a reduction of domain area, 
figure [5] shows that the narrower the domain, the larger the drag reduction ratio. Within the uncertainty the 
behavior is almost linear and is most probably due to a combination of a virtually further reduced domain 
size in viscous units in case of high DR, bringing along an underestimation of C/, as well as to a more 
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FIG. 6. Percentage drag reduction lOOi? versus oscillation period T + for the oscillating wall, at A + = 12. Blu open 
circles: Re T — 200; red open squares: Re T = 1000. Error bars at 95% confidence level. Reference DNS data: blu 
solid line and filled circles, Re T = 200 from2£; filled green squares, Re T = 400 from2fi; open green squares: Re T = 500 
froro^; closed red square: Re T = 1000 from22. 



ordered structure of turbulence, prone to be affected by the control. The latter is both a drawback and 
a potential strength of the method, which could give a more comprehensible picture of the drag reduction 
mechanism, by discarding part of the turbulent structure that do not directly participate in the autonomous 
regeneration cycle and offering a more ordered array of elemental turbulent structures. Nevertheless, it is 
important to point out that the reported case of figure [5] is in the region of the parameter space that is 
mostly affected by this effect. When traveling waves are employed and, in general, when the drag reduction 
is below 0.3, this effect is not observed. 



III. DRAG REDUCTION RATIO 

Figure [6] represents the drag reduction ratio R for a wall oscillating with a non-dimensional period T+. 
Here and throughout the rest of the paper our new direct numerical simulation are presented as big empty 
colored symbols. Color and shape refer to different values of the Reynolds number: blue circles at Re T = 200, 
red squares at Re T = 1000 and green triangles at Re T = 2000. Error bars have been added to the presently 
computed data, for the 95% of confidence level. (The error bars are not seen when their size is within the 
marker). Literature data at different Re are drawn with other colors and specified in the legend. Our data 
are superimposed to well-established, low-i?e, DNS datasets represented as filled circles, in particular those 
at Re T = 200 by Quadrio and Ricco^i and those at Re T — 400 by Ricco and Quadrio^P. for the oscillating 
wall (k x = 0), and those at Re T = 200 by Quadrio et al3& for the traveling waves. When available, higher- Re 
simulations are reported: here those at Re T = 500 and Re T = 1000 by Touber and Leschziner— . 

The already available DNS data span a significant part of the parameter space only at the low Reynolds 
number of Re T = 200, while the more recently computed data points at higher values of Reynolds number 
sample the parameter space more sparsely and focus only on the optimal period T + m 125 or on its 
neighborhood. 

We observe that the values of R obtained with the present simulations at Re T = 200 agree very well 
with those obtained with full-channel simulations in Quadrio and Ricco^I at the same Reynolds number. 
The agreement is perfect at values of T + much larger and much smaller than the optimum value. Near 
the optimum, though the position of the maximum is well-captured, a slight ovcrcstimation of the drag 
reduction occurs, together with an increase of the uncertainty. This overestimation decreases to only 1% 
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FIG. 7. Percentage drag reduction ratio 100i? versus oscillating frequency lj + for the streamwise-travelling waves, 
at A + = 12 and A + = 1250. Blu open circles: Re T = 200; red open squares: Re T — 1000; green open triangles: 
Re T = 2000. Error bars at 95% confidence level. Reference DNS data: blu solid line and filled circles, Re T — 200 
from2£. 



for the traveling waves. In the followings reference full-channel DNS are reported to confirm the very good 
agreement with present simulations at low Re that will not be underlined again. It is interesting to point 
out the agreement with the point taken from Touber and Leschziner— at Re T = 1000, though it lies on 
slightly bigger value of DR. The curious overestimation of Touber's data is even more evident at Re T = 500: 
their point at T + = 200 has an higher DR than the point obtained via full-channel DNS at Re T = 400. 

When the Reynolds is increased to Re T = 1000 our dataset shows that the maximum drag reduction ratio 
consequently reduces from 0.39 to 0.29, and occurs at T + = 90, thus the data are consistent with a shift of 
the optimal period towards lesser values. The shift is a result of the unevenly sensitivity to Re throughout 
the whole range of T+. In particular, as suggested by Ricco and Quadrio^, the decrease appears to be much 
less significant at low values of T + , while it is amplified at high T + . As a result a similar DR is achieved at 
T+ = 50 and T+ = 100. 

This scenario of heterogeneous dependence of R upon Re is fully confirmed if we enlarge to the k~£-lo + 
plane of the traveling waves at Re T up to 2000. Figure [7] plots lOOi? for the horizontal slice numbered (2) 
in figure [3] i.e. as a function of the forcing frequency w + for waves with constant wavelength A + = 1250. 
The maximum drag reduction ratio achieved at this Re is very high and equal to 0.49, while the reference 
data provide R = 0.48, with a difference clearly within the confidence interval. 

As the Reynolds number is increased to Re T = 1000, K ma x drops to 0.37, with a total loss of 0.12. A 
similarly significant change can be observed in the "valley" of R, where at low Re the waves produce a 
negative R, i.e. a drag increase. Although the drag increase almost disappears, the minimum of R is still 
present and seems to widen, embracing a larger range of w + . Again, the data suggest that the optimal value 
of the control parameters might change with Re. At this wavelength, the optimal point is the standing 
wave, bj + = 0, at higher Re. As in the previously discussed case of the oscillating wall, the most important 
observation that can be drawn from the figure is that the influence of Re is heterogeneous over the parameter 
space, being stronger approaching the region of DI. Outside of this region, and in particular for \uj + \ > 0.15, 
R remains nearly unchanged. 

The dataset at Re T = 2000 consists of four points only that indeed confirm the general trend discussed 
above: the DI valley is widening and influencing the neighborhoods. It pushes high-DR regions on its right 
toward higher frequencies and on its left toward lower frequencies. The highest DR on the slice is R = 0.29, 
at uj = -0.05. 

In figure [51 the interest is focused on the low-i?e ridge of maximum R, i.e. the locus of points with 
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FIG. 8. Percentage drag reduction lOOii versus the oscillating frequency u) + for the streamwise-travelling waves, at 
A + — 12 and along the ridge of the low- Tie maximum, defined as c + = 2.5. Symbols as in figure]?] 
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FIG. 9. Percentage drag reduction 100R versus the wavenumber k+ for the streamwise-travelling waves, at A + = 12 
and uj + = 0.012. Symbols as in figure [Jj 



maximum R that approximately extends along the line c + 2.5 (line (3) in figure [3]). The low-frequency 
low-wavenumber ('southern') region of the ridge is particularly affected by a change in Re while the right 
high-frequency part is almost insensitive: here we encountered at Re T = 1000 the maximum R of 0.42. The 
DR ridge seems to be moving upwards or, in other words, its 'southern' 1ow-k+ tail is receding. 

This scenario is fully confirmed by the two further constant-a; + sections of the k — uj plane, namely at 
oj + = 0.012 and uj + = —0.2, shown in figure l9l and fTOl respectively. From the first the receding behavior 
of the low-wavenumber part of the ridge is extremely evident: at k + = 0.005 the DR almost halves by 
increasing Re r to 1000. In the latter the trend of R towards higher wavenumbcrs is milder but still present, 
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FIG. 10. Percentage drag reduction 100i? versus the wavenumber for the streamwise-travelling waves, at A + = 12 
and lo + = -0.2. 
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FIG. 11. In the filled circles the value of the Reynolds sensitivity exponent 7, supposing R ~ Rel superimposed to 
the drag reduction map aftes2£. 



though shaded by the large confidence interval. 



A. Dependence upon Re 



To quantify the i?e-sensitivity, now and in the following we will adopt the usual supposition that in the 
present Reynolds number range R is proportional to Rel . This common choice is rather arbitrary and based 
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FIG. 12. Input Power Pi n versus wave frequency uj + at various Re at A + — 12 and A + = 1250. Open symbols and 
colors as in figure [3 Lines are the laminar solution at Re T = 200 (blue), Re T — 1000 (red) and Re T = 2000 (green). 
The magenta line is the phase lag between the spanwise wall velocity w w and its derivative dw/dy from the laminar 
solution. 

on low-i?e observations but it is useful as a benchmark in this range of Re, even though other choices have 
been proposed or can be thought, an example^ is R = a + Re@, where (3 is eventually f3(Re). 

It is clear from the present study that the sensitivity coefficient 7 is not unique. In the case of oscillating 
wall (k x = 0), the empirical coefficient 7 w —0.20 is confirmed for the maximum drag reduction R,„ al . The 
slightly higher value of 7 = —0.25 is achieved at constant T + = 100 because the optimal period is shortening 
with Re. On the other hand, drag reduction ratios at short periods of forcing, for example corresponding 
to u + = 0.2 are less sensitive to a change in Re, leading to 7 « —0.08. 

If we extend to the traveling waves, the scenario is unexpectedly rich. The drag reduction ridge is reduced 
in size and looses its low-frequency low-wavenumber tail, where we locally measure a high sensitivity of 
7 = —0.25. The trend suggests that this part of the drag reduction ridge is a low-i?e effect, that is deemed 
to disappear quickly with Re. The high-wavenumber high-frequency part is instead less sensitive ans shows 
a 7 = -0.08 at (lu+ = 0.04, k+ = 0.015) which reduces to 7 = -0.07 at (uj+ = 0.08, n+ = 0.035). This 
observation is of great importance since the upper part of the ridge is still a region of high drag reduction 
ratios, ranging between 0.4 and 0.5, suggesting the possibility to preserve significant value of R at high Re. 
Surprising is also the behavior of the DI region: its intensity is reduced while its extension seems to widen 
and influence the neighborhood. At Re T = 1000 the drag increase almost disappears or turns to very little 
DR. In the latter case a sensitivity exponent can not even be defined within the framework of the R ~ ReJ 
law and around (w + = 0.06, re+ = 0.005) goes from 3 to 00. Finally we observe that at high frequencies 
(|cj + | = 0.2) both forward- and backward-traveling waves are not strongly affected by a change of Re: 7 is 
-0.08 at low wavenumbers and decreases further to -0.05 at k+ = 0.04. 

To summarize, the i?e-effect is proved to be lower at higher wavenumbers and at increasing distances from 
the DI region, whose influence increases with Re, though its intensity reduces. This implies that on the left 
of the DI valley the more insensitive regions move towards lesser frequencies, even negative when necessary, 
and higher wavenumbers, while the opposite occurs to the other side. As a consequence, the optimal set of 
control parameters is not constant with Re. 

IV. INPUT POWER AND POWER BUDGET 

From figure [12] one sees first that the laminar GSL solution predicts the input power rather well, but 
regions of the parameter space exist where prediction is not perfect. In particular when the equivalent 
period, defined by Quadrio et &&■ as T + = A+/(M+ — c+) is about ten. U w is the convective speed of 
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FIG. 13. Map of the change AS = S(Re T = 1000) - S(Re T = 200) in the net power budget S: in the circle the value 
AS. Arrows points the value of S at Re T = 1000. Blue circles: negative S; Red circles: positive S. 

the skin-friction information, as defined in Quadrio & Luchini^i . In this region, the traveling-wave is 
seen as standing by the near-wall turbulent structures with which shares the same velocity, leading to a 
lock-in effect that leads to drag increase and deviations from the laminar solution of the generalized Stokes 
layer. The second, most important observation is that Pm/Po, that is the percentage of the pumping power 
required to run the flow, slightly decreases for increasing Re. Another way of putting it is that Pi n is 
constant when expressed in inner units. This was suggested by Ricco & Quadrio^, can be verified through 
a nondimcnsionalization of the laminar expression of Pi n and is fully confirmed here. Whenever the laminar 
solution holds, the input power is seen to decrease as Pe~ 0136 . This opens up interesting perspective as far 
as the net energy saving S is concerned: in several parts of the control parameter space Pj n /Po decreases 
faster R, meaning that S could even increase with Re. Evidence is given in figure H31 where the change 
AS in net power budget when Re T is increased from 200 to 1000 is represented. At these A + most of the 
parameter space at Re T = 1000 still has S < (blue circles), though rapidly increasing and sometimes even 
approaching positive values. For control parameters lying on the low-wavenumber tail of the DR ridge S 
is positive (red circle) but showed a decrease of 7%. Here R decreases faster than Pj n /Prj- O n the upper 
part of the ridge the trend is opposite, a positive and increasing S is achieved. This result is encouraging 
especially if we consider that much lower A + should be considered when interest is focused on the net power 
budget. 



V. CONCLUDING DISCUSSION 

The main messages conveyed by the results presented so far is that the general picture, as it stands today, 
of how the drag reduction ability of techniques based on spanwisc forcing decreases as Re increases might 
need to be revised. Our data indicate first of all that that drag reduction decrease with Re at a rate which 
depends strongly on the region of the parameters space. Therefore, the optimal set of control parameters 
changes with Re and new regions of the parameter space may lead to interesting values of drag reduction 
and net power saving, meaning that a meaningful parametric study has to span a large portion of parameters 
space and not only focus on the neighborhood of the low-Pe optimum. 

A possible scenario is shown in figure 1151 The optimal control parameters shift towards higher wavenum- 
bers and away from the drag increase region, whose intensity strongly reduces with Re, while the portion 
of parameter space affected by its influence seems to increase. The Pe-effect appears to both reduce drag 
reduction and widen the DI region of influence, suggesting that two still unknown different mechanisms 
exists through which the Pe-effect acts. 
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FIG. 14. Relative change in R as function of Re T for traveling waves with wavenumber kJ" = 0.005. On the left 
lo + = —0.2, on the right cj + = 0.005. Both results with the reduced channel and bespoke fullchannel simulations are 
represented. 



The power law R ~ Re" traditionally used to quantify the .Re-effect may be too simple a law to clearly 
depict the complex scenario and clearly fails in cases where the drag increase (RjO) turns to drag reduction 
(RiO). 

We are well aware of the limitations implicit in the present study, which in particular has its weakest 
spot in the rather limited size of the computational domain. Although first-order effects can be ruled out, 
since the domain is always many many times larger than the MFU, the present results should be considered 
as suggestive only, though predictions quantitatively agree with simulations of Touber & Leschziner— . 
Moreover, they are in line with other full-DNS data recently obtained by Hurst, Chung & LockerbyH. 

On the other hand, if our simplified framework well represents the drag reduction, it would imply that 
large scale motions only marginally participate in the mechanisms behind the i?e-cffcct which indeed would 
be due to near-wall dynamics. This scenario is even more interesting considered that dynamical importance 
of linear phenomena in the near-wall region, giving hope for simplified representations of drag reducing 
mechanisms and confirming the studies of Duque et alJ^. 
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